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Abstract. 

Above a critical velocity, the dominant mechanism of energy transfer between 
a moving object and a dilute Bose-Einstein condensate is vortex formation. In this 
paper, we discuss the critical velocity for vortex formation and the link between 
vortex shedding and drag in both homogeneous and inhomogcneous condensates. 
We find that at supersonic velocities sound radiation also contributes significantly 
to the drag force. 



PACS numbers: 03.75.Fi, 67.40.Vs, 67.57.De 
1. Introduction 

One enticing consequence of the discovery of Bose-Einstein condensation (BEC) in 
dilute alkali vapours Q is the potential for refining our understanding of quantum 
fluids. In particular, the dilute Bose gas provides a near-ideal testing ground for 
elucidating the role of vortices in the onset of dissipation in superfiuids. Recent 
experiments have demonstrated the formation of quantized vortices by rotational 
excitation of one Q and two-component condensates || , in analogy with the famous 
'rotating bucket' experiments in liquid helium In addition, by moving a far-off 
resonant laser beam through a condensate, Raman et al. Q measure a heating rate 
which suggests a critical velocity for dissipation characteristic of vortex shedding. 

The attractive feature of experiments on dilute Bose gases is that the weakly- 
interacting limit permits quantitative comparison between theory and experiment. 
The theoretical description is based on the Gross-Pitaevskii equation, a form of non- 
linear Schrodinger equation (NLSE) ||. In the NLSE model, the critical velocity for 
vortex shedding by a cylindrical object was found to be, v c ~ 0.4c, where c is the 
speed of sound 0. In a trapped (inhomogeneous) condensate the critical velocity is 
lower due to the reduction in density, and hence sound speed, towards the edge of the 
condensate [[| ||. 

In this paper, we study vortex shedding due to the motion of an object through 
homogeneous and trapped Bose-Einstein condensates using the NLSE model. The 
hydrodynamical properties of a NLSE fluid are reviewed in Sec. In Sec. ^ and 
we discuss the critical velocity for vortex nucleation and the link between vortex 
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formation and drag in homogeneous condensates. In Sec. || we consider the properties 
of trapped condensates and highlight the differences from the homogeneous case. 

2. Quantum fluid mechanics 

At low temperatures and low densities, atoms interact by elastic s-wave scattering, 
and collisions can be parameterised by a single variable, the scattering length, a. For 
atoms of mass m, the wavefunction of the condensate, ip(r, t), is given by the solution 
of the time-dependent Schrodinger equation: 
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where the wavefunction is normalised to the number of atoms, AT, the coefficient of 
the non-linear term, g = ATrfi 2 a/m, describes the interactions within the fluid, and 
V(r,t) represents external potentials arising from the trap and any moving obstacle. 



2.1. Fluid equations 

The link between the nonlinear Schrodinger equation (|l|) and the equivalent equations 
of fluid mechanics is well known ^]. However some points concerning condensate 
flow near penetrable objects are less familiar. In this section we gather some of the 
key concepts and equations that figure in the discussion of our simulations. 

Classical (isentropic) fluid mechanics is based on two coupled differential 
equations: one describing the transport of mass, the other the transport of momentum 
[|l0| . The relevant quantum variables can be constructed from the wavefunction: the 
mass density p and momentum current density J are defined as, 

p ee mi/j*^ and J k = (h/2i){ip*d k iP - ^d k ^*) , (2) 

where the index k denotes the vector component. The fluid velocity is defined 
by Vk = Jk/p, or equivalently in terms of the phase, (f>, of the wavefunction, 
Vk = (fi/m)dk<t>. Clearly, the velocity field is a potential flow, however it is also 
compressible and furthermore can support circulation (vorticity) as will be seen. 

The conservation of mass (probability), i.e., the continuity equation, follows from 
the definition of p and equation (|l|) 

d t p + d k J k = 0- (3) 

The conservation of momentum equation may be found by considering the rate of 
change of the momentum current density, 

d t Jk + djT jk + pd k (V/m) = , (4) 

where the momentum flux density tensor takes the form 0, [llj , 

T jk - ^-(djpd k 1> - rdjd k i> + c.c.) + l^fel^l 4 ■ (5) 
This can be rewritten as, 

Tjk = pVjVk - Ojk , (6) 
where the stress tensor o~jk is given by, 

°j* = -¥]k9{p/m) 2 + (h/2m) 2 pd 3 d k hip ■ (7) 
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2.2. Pressure, sound and drag 

The form of equations (||), (|j), and (||) is identical to those for classical fluid flow 
(|To|l , the difference emerges from the nature of the stress, equation (0). A classical 
ideal fluid is characterised by ajk = 0, for all j, k. In a viscous fluid, the shear stress 
( a jk , j 7^ k) is produced by velocity gradients between neighbouring streams such that, 
ajk — r](djVk + dkVj), where r\ is the coefficient of viscosity. This creates a frictional 
force which gives rise to energy loss. In a pure dilute Bose-Einstein condensate there 
is no frictional viscosity, but a shear stress arises from density gradients, the second 
term in equation (Q). This property gives rise to the possibility of vortex formation 
and drag without viscosity. 

The pressure (normal stress, —o~jk,j — k) within the quantum fluid takes the 
simple form, 

p = \g{p/mf - (h/2m) 2 pV 2 In p . (8) 

The second term, called the quantum pressure, is weak in homogeneous regions of the 
fluid, that is, far from obstacles or boundaries, vortex lines or shocks. The essential 
difference between interacting and noninteracting (ideal) fluids is the existence of 
interaction pressure which supports sound propagation. In the bulk of the fluid, where 
the quantum pressure is negligible, the speed of sound || H, 

c= Vdp/dp= (gp/m 2 )i . (9) 

The force on an obstacle moving through a condensate can be calculated from 
the rate of momentum transfer to the fluid. By integrating Eq. (0) , one finds that the 
fc-th component of the force, 

F k = d t I dfijfc = - [ dS nj T jk - I dfl pd k {V/m) , (10) 
Jn Js Jn 

where S is the surface of the object or control surface within the fluid jn]], fl is the 
volume enclosed by S, rij is the j-component of the normal vector to 5, and dS is 
a surface element. The second term on the right-hand side can be likened to the 
buoyancy of the fluid. In the case of homogeneous flow past an impenetrable object 
(Sec. |] and ^), the wavefunction vanishes on the object surface and the potential 
is uniform elsewhere, therefore, only the first term contributes. Conversely, for a 
penetrable object in a trapped condensate (Sec. ||), Q may be chosen to encompass 
the entire fluid, and the first term is negligible compared to the second. 



2. 3. Quantisation of circulation 

The quantum Euler equation follows from combining the equations describing the 
conservation of mass and momentum, (|^) and (|]), along with the identity, 

p^djipd.dk inp] = 2d k \p-^d j a j p^] , (n) 

allowing the momentum equation to be written as, 

d t Vk + vjdjVk + d k [gp/m 2 - (h 2 /2m)p-^d 3 djpi + V/rn] = . (12) 

The conservation of energy (Bernoulli equation) then follows as an integral of Euler's 
equation, or more directly from the real part of equation ([!]): 

hd t (f> + \mv 2 + gp/m - (h 2 /2m)p-^ V 2 ^ +V = 0. (13) 
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Perhaps the most significant quantum effect on the mechanics of the fluid is the 
quantisation of angular momentum. The circulation is given by, 

r = <j)dr ■ v = (h/m)2ns s = 0,1,2,... (14) 

where the closed contour joins fluid particles. The conservation of angular momentum 
(Kelvin's theorem), follows from Euler's equation jl^ ) and states that the circulation 
around a closed 'fluid' contour does not change in time. This means that within the 
fluid, vortex lines must created in pairs which emerge from a point. The exception 
is at boundaries, where the wavefunction is clamped to zero and no closed fluid loop 
can be drawn, e.g., at the surface of an impenetrable object jio) or from the edge of 
a trapped condensate |ll| . 



2.4. Units 



For a homogeneous fluid flow, where the external potential is due to the obstacle 
only, it is convenient to rescale length and velocity in terms of the healing length 
£ = fij \jmnQg and the asymptotic speed of sound c ~ y/ n^g/m, respectively. In this 
case, equation (0) becomes 



idtij(r,t) = 



V(r') + \ij(r,t)\ 2 j(r,t) , 



(15) 



where ip — ■0/ v / «o & n d no is the number density far from the object. The force per 
unit length is measured in units of Ti^ri^g /m. Unless otherwise stated we use these 
units throughout. For steady flow, in which v and n are independent of time and 
(j) = —j-it, where [i is the chemical potential, the Bernoulli equation takes the form 

- V + \v 2 = constant . (16) 



3. The critical velocity 

The critical velocity for the breakdown of superfluidity is often expressed in terms of 
the Landau condition v c = {e/p) m - m , where e and p are the energy and momentum 
of elementary excitations in the fluid, and v c is the flow velocity in the fluid bulk . 
In the dilute Bose gas, the long wavelength elementary excitations are sound waves 
and the Landau criterion predicts that v c = c. However, for flow past an object, the 
local velocity near the obstacle, v, can become supersonic even when flow velocity, U, 
is subsonic. Consequently, the critical flow velocity, v c , where laminar flow becomes 
unstable occurs at a fraction of the sound speed. An estimate of v c may be found 
following the argument suggested by Frisch et al. . For an incompressible flow past 



a solid object, Bernoulli's equation (16) (neglecting the quantum pressure) has the 
simple form, 

n(v) + \v 2 = 1 + \U 2 , (17) 

where U is the background flow velocity. The maximum velocity which occurs at 
the equator of the object is v = |J7 for a sphere (or v — 2U for a cylinder), 
therefore, ^(| — 1)U 2 = 1 — n{v). The critical velocity is reached when the maximum 
speed, v = is equal to the 'bulk' sound speed, c = y / n(w), which gives 

v c = U = -\/8/23 « 0.59 (or y / 2/ll w 0.44 for a cylinder). However, for a compressible 
fluid, the equatorial velocity is slightly larger due to pressure effects. The first-order 
correction gives v = |[7 + ^U 3 which reduces the critical velocity to v c « 0.53. 
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Figure 1. Laminar flow past a sphere with radius R = 50 for two flow velocities 
U = 0.2 (left) and U = 0.53285 (right). The three curves show, (a) the 'bulk' 
sound speed c = y/n = \ip\, (b) the velocity, v, and (c) the quantum pressure, 
V 2 i/ , /2i/', as a function of position. Note that, at the critical velocity u = 0.53285 
(right), the fluid velocity is equal to the sound speed where the quantum pressure 
is exactly zero. 

The exact value of v c may be found by solving the uniform flow equation [l^] , 

id t 4>'(r',t) = -§V' 2 + 7(rO + |V>'(r / ,i)| 2 + h>- V]tp'(r',t) , (18) 

where ij)'(r',t) — ij)(r,t) is the wavefunction in the fluid rest frame written in terms 
of the object frame coordinates, r' = r — vt. Stationary solutions of the form, 
ip'(r',t) = <fi(r')e~ Hlt are found to exist only for v < v c , where v c is the critical 
velocity for vortex formation. 

To illustrate the behaviour of the exact solutions near the critical velocity, we solve 
equation (fsl) in 3D for an impenetrable sphere with radius R — 50. The wavefunction, 
velocity and quantum pressure term near the object are shown in Fig. |^. Note that 
these parameters are related via the Bernoulli equation ( |l6| ) . The intersection between 
the velocity v and wavefunction amplitude, \ip\, curves defines the position where the 
velocity is equal to the 'bulk' sound speed (||). Note that close to the object the 
effective sound speed is increased due to the quantum pressure term, (||), therefore 
even though the density is low the flow is not 'supersonic'. The critical velocity is 
reached when flow velocity exceeds the speed of sound in the bulk of the fluid, i.e., 
when the intersection between the velocity and wavefunction curves moves into the 
region where the quantum pressure term is zero, see Fig. |l|(right). 

The complication for trapped condensates is that the speed of sound and hence 
the critical velocity depend upon position. In addition, the object potential is typically 
penetrable and non-uniform. We return to these topics in Sec. |5|. 

4. Vortex shedding, drag, and dissipation 

For flow faster than the critical velocity, v > v c , vortices are emitted approximately 
periodically. A typical vortex stream pattern for flow past an impenetrable cylinder is 
shown in Fig. ||. The background flow is from right to left and the vortex - anti- vortex 
pairs produce a flow pattern which opposes the background. Consequently, the vortex 
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Figure 2. A shaded contour plot showing the condensate density for 
homogeneous flow past a cylinder with radius R = 3. The dark regions indicate 
lower density, i.e., the position of the object and vortex lines. 

trail separates the main flow from an almost stationary wake. The momentum loss 
from the fluid is transferred to the object creating a drag force. The contribution 
of vortex shedding to the drag force can be estimated by considering the momentum 
transfer due to vortex emission, i.e., 

F v = f vPv , (19) 

where f v is the vortex shedding frequency and p v is the momentum of a vortex pair 
as it is created at the equatorial plane. Small fluctuations in the vortex shedding 
frequency occur because as the vortices move downstream they interact with each 
other creating fluctuations in the flow pattern around the object. The drag force is 
taken to be the time-average over many vortex emission cycles. 

In addition to vortex shedding, drag may also arise due to sound waves. The 
time-average drag, evaluated using equation (|l0|), as a function of velocity for an 
impenetrable cylinder is shown in Fig. |^. The drag is zero up to the critical 
velocity, then increases approximately quadratically with v (?], [yj. Also shown is 
the contribution to the drag force from vortex shedding alone, i.e., equation (|l9|). 
This comparison illustrates that for v < c the drag is produced by vortex shedding, 
whereas for v > c, an increasingly significant contribution arises from sound waves. 
For v > c, the reflected matter waves create a standing wave pattern in front of the 
object as shown in Fig. [I| 

Using similar arguments for energy loss of the flow. Consider now the condensate 
at rest and the obstacle moving with speed v. A drag force leads to energy transfer 
to the condensate. The energy transfer rate is given by 

AE , . 

= ^drag-U ■ (20) 

Eqs. (|l^) and (^0|) make the important link between vortex shedding, drag and energy 
dissipation. 

5. Motion in a trapped condensate 

The inhomogeneous density profile and finite size of trapped condensates means that a 
steady, uniform flow is difficult to achieve. The MIT experiment [^| partially overcame 
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Figure 3. The time-averaged drag force as a function of velocity for an 

impenetrable cylinder with radius R = 3. The open circles indicates the 

contribution due to vortex shedding. The error bars indicate the residual 
fluctuations in the time-averaged drag. 




Figure 4. Time-averaged density for homogeneous flow past a cylinder (radius 
R = 3) for flow speeds of 0.9 (left) and 1.2 (right) times the speed of sound. For 
supersonic flow a standing wave pattern appears in front of the object. 
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Figure 5. The drag force as a function of velocity for a laser beam moving back 
and forth through an inhomogeneous condensate. The triangles, squares and 
circles corresponds to 2D simulations with different oscillation amplitudes. The 
diamonds correspond to a 3D simulation. In 2D the force is effectively that on a 
disk of depth 2.5 healing lengths, whereas in 3D the condensate depth corresponds 
to a Thomas Fermi radius of 6.4 healing lengths. 



this problem by sweeping the object back and forth at constant velocity within the 
central region of the condensate, where the density is approximately uniform. In this 
case, the object moves through its own low-density wake, and consequently the drag 
law is different from the uniform flow case discussed above. 

In Fig. | we show the time-averaged drag on a laser beam oscillating in a trapped 
condensate. Important differences with the uniform flow case (Fig. ^) arise at high 
and low velocities. The drag force tends to saturate at higher velocities as the object 
expels fluid from the region of oscillation and the pressure drops. 

Fig. ^ also highlights important differences between two and three dimensions. 
Two dimensions corresponds to the limit of a 'cylindrical' condensate, where the 
density is uniform parallel to the axis of the object (which wc dehnc as the z-axis). 
However, in realistic three-dimensional situations the density is inhomogeneous along 
z, leading to a variation of the speed of sound which vanishes at the condensate edge. 
This results in a lower critical velocity in 3D than in 2D, as apparent in Fig. [| In the 
Thomas-Fermi limit (go? 3> 1) the density profile is parabolic, therefore the average 
sound speed along z is a factor of 7r/4 smaller than that at the centre. The remaining 
reduction arises from the fact that even at very low velocities, vortices tend to be 
formed where the laser beam intersects the lower density fluid at the condensate edge. 
The lower densities arising in 3D also leads to enhanced sound emission, and hence 
enhanced drag at high velocities. 

Below the critical velocity, dissipation due to sound emission occurs at the motion 
extrema, where the object accelerates. This is illustrated in Fig. || which shows a cross 
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Figure 6. Cross-section through a 2D trapped condensate pierced by a Gaussian 
laser beam 'object' (also shown), (a) If the object is moving at constant velocity, 
the fluid density is symmetric around the object, and the net drag is zero, (b) 
If the object accelerates the fluid cannot adjust sufficiently rapidly, leading to an 
asymmetric fluid distribution around the object which produces drag. The fluid 
relaxes by the emission of a sound wave. The a;-axis is in units of the laser beam 
waist, while for the y-axis the units are arbitrary. 

section through a 2D condensate cut by a moving laser beam. For constant motion, 
Fig. ^](a), the fluid is distributed symmetrically around the object and the drag, which 
is given by an overlap integral between the condensate density and the gradient of the 
object potential (i.e., the second term in equation (fiol)) is zero. When the object 
accelerates, Fig. |g(b), the fluid fails to respond rapidly enough to the abrupt change 
in velocity, and the asymmetry in the overlap between the fluid and the objects leads to 
a resistance force analogous to dynamic buoyancy. The system relaxes to the uniform 
flow case by the emission of a sound wave. This corresponds to the 'phonon heating' 
process discussed in ||. 

Although the NLSE model describes dissipation in the sense of energy transfer 
between a moving object and the fluid, it does not say anything about how that energy 
(mostly stored within the vortex core) may be subsequently converted into heat. A 
complete description should include coupling of the condensate to a thermal cloud, 
and would describe the damping of phonon and vortex modes. Recent work on the 
non-equilibrium dynamics of the condensate and non-condensate predict a depletion 
of the condensate fraction Jl(| as observed in the MIT experiment. 

6. Summary 

The motion of an object through a dilute Bose-Einstein condensate provides an ideal 
system to study the fundamental problem of the onset of dissipation in superfluids. In 
this paper we have explained the role of vortex shedding and sound emission in energy 
transfer between the object and the condensate. No energy transfer is observed under 
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the condition of uniform, steady flow at speeds below a critical velocity. However, if 
the object accelerates there is a small dissipativc effect, even below the critical velocity, 
due to sound emission. The critical velocity is reached when the local velocity in the 
bulk of the fluid (i.e., where the quantum pressure is zero) exceeds the speed of sound. 
Above the critical velocity vortices are emitted leading to a drag force and energy 
transfer to the fluid. Vortex shedding dominates the energy transfer for intermediate 
velocities, while sound emission becomes increasingly important for supersonic motion. 
We highlight some important differences between homogeneous and inhomogeneous 
trapped condensates. In particular, that the critical velocity is substantially reduced 
when the object intersects regions of lower density at the condensate edge and that 
the trap inhomogeneity gives rise to an additional term in the drag force analogous to 
a buoyancy. 
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